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We reproduce apparently complex cellular automaton behaviour with simple partial differential 
equations as developed in Our PDE model easily explains behaviour observed in selected sce- 
narios of the cellular automaton wargame ISAAC without resorting to anthropomorphisation of 
autonomous 'agents'. The insinuation that agents have a reasoning and planning ability is replaced 
with a deterministic numerical approximation which encapsulates basic motivational factors and 
demonstrates a variety of spatial behaviours approximating the mean behaviour of the ISAAC sce- 
narios. All scenarios presented here highlight the dangers associated with attributing intelligent 
reasoning to behaviour shown, when this can be explained quite simply through the effects of the 
terms in our equations. A continuum of forces is able to behave in a manner similar to a collection of 
individual autonomous agents, and shows decentralised self-organisation and adaptation of tactics 
to suit a variety of combat situations. 
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I. INTRODUCTION 

Lanchester modelled force dynamics as a set of ordi- 
nary differential equations that have greatly influenced 
military decision making for many years and permeate 
military thought and analysis to this day. However, the 
absence of any spatial representation has been a con- 
stant source of criticism for using any form of the Lanch- 
ester Equations. Extensions have been made by Pro- 
topopescue et al. who treated the spatial aspect using a 
simple advection-diffusion representation as is typical in 
predator-prey type modelling 0, H, 0] ■ Despite using spe- 
cific initial density profiles and short simulation times to 
prevent excessive and unrealistic diffusion of the forces, 
this approach was able to represent a suite of realistic 
behaviour - frontal attack, envelopment, infiltration and 
turning manoeuvres. However, replacing the constant ve- 
locity with a spatially dependent velocity field resulted 
in unacceptable numerical losses, restricting the velocity 
field to a temporally dependent one. 

In Pi a set of partial differential equations for a more 
spatially realistic approach to combat modelling was pro- 
posed. Through the implementation of biological aggre- 
gation models, a series of basic scenarios demonstrat- 
ing cohesive soldier movement throughout the domain 
regardless of losses incurred through firing effects were 
explored. It was also shown that a scenario can be 
viewed as a series of sub-battles characterised by peri- 
ods of constant density loss which was of the same form 
as in 0, H, 01 • This was the intended method of imple- 
mentation of the LEs as suggested by Lanchester himself. 
These PDEs were developed with the intention of provid- 
ing an alternate method for investigating the dynamics 
seen in cellular automaton-based wargames. 
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These agent-based or cellular-automaton (CA) models 
have received much attention in many disciplines, partic- 
ularly in defence related research. Models such as Ein- 
stein [5*1 , ISAAC , Map Aware Non-uniform Automata 
(MANA) demonstrate a range of behaviour which ap- 
pears to hint at some form of underlying structure. Each 
individual troop is modelled via a rule set relating to 
quantifiable capabilities such as fire-power, communica- 
tions and also intangibles such as morale or desire to 
remain close to friendly forces. These rules encode the 
nonlinearities necessary for a more realistic description 
of warfare. These nonlinearities need to be understood 
in order to develop specialised tactics based on current 
capability, or enhance the procurement of future capa- 
bility. Ilachinski [8], who has been instrumental in the 
development of ISAAC, stresses the need for research 
into nonlinear continuous dynamics, exploitation of anal- 
ogous biological models and phase-space reconstruction 
techniques. Lauren compares MANA simulation results 
with fluid dynamic concepts or transition between lami- 
nar and turbulent states and maintenance of force profiles 
to viscosity [9,] . 



Continuous models can be more transparent in terms 
of how parameter changes affect outcomes and are thus 
more understandable. In this paper, the PDE model is 
expanded to refiect the underlying assumptions of the 
ISAAC model to provide an alternate analysis of the 
demonstrated behaviour. By effectively removing the 
inherent stochasticity of ISAAC, we can determine the 
mean behaviour of a scenario and more easily under- 
stand the effects of parameter changes. Introducing spa- 
tial asymmetry through the initial density profile position 
forms a type of controlled stochasticity approximation, 
providing an explanation of the observed ISAAC dynam- 
ics. 
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A. PDE model 

As developed in [ij, we begin with the following 
integro-differential equations in two dimensions. Each 
equation is separated into three parts, fdiff containing 
the diffusion terms, fvei containing the velocity terms 
and freact Containing the interaction terms. 

Ou 

— ■ (Du(m)Vw) + u{ku*v) + duV + 
ot ^ ^ ^ ^ J ' 

fdifF freact 

V • {u(CuU + Aa{Ka *U)- ArU{Kr * u))} (f ) 



dv 

— = V • (Dv(w)Vu) +v{k^*u) + d^u + 
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fdiff freact 

V • {v{C^V + Aa{Ka * V) - ArV{Kr * v))} (2) 



Attraction and repulsion operate over different spatial 
domains, and respectively with > r^, and a spa- 
tially dependent vector field C is used. The form of the 

kernel in fr, 
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Following the work of Boswell [13] and Hundsdorfer 
[llj l , the Method of Lines scheme is implemented where 
the spatial derivatives are discretised, giving a large sys- 
tem of Ordinary Differential Equations. Flux limiters are 
employed at this stage to ensure positivity and conser- 
vation of mass. An explicit Runge-Kutta method is then 
used for the time integration, again with constraints in 
place to ensure positivity. 



II. ISAAC 

In this section the cellular automata model ISAAC 
and a selection of scenarios its authors have published to 
demonstrate its range of capabilities are described. We 
then show how our simulation with slight modifications 
can reproduce similar behaviours. Thus we show that a 
PDE model can perform the combat modelling tasks of 
a state-of-the-art CA model while being easier to under- 
stand and analyse. This can have great advantages when 
seemingly unexplainable or novel behaviour is seen in CA 
modelling results that may otherwise be attributed to a 
form of intelligence. 

ISAAC, Irreducible Semi- Autonomous Adaptive Com- 
bat, is a multi- agent based simulation in the style of a 
cellular automaton. A CA is essentially a lattice where 
information located on the nodes propagates from its po- 
sition at each time step based on a set of defined rules. 
Ilachinski asserts ISAAC differs from a traditional CA 
in that agents rather than information move throughout 
the lattice and that rule sets may adapt over time. In the 



comparison ISAAC scenarios, however, rule sets remain 
fixed. 

In ISAAC, each agent represents a simplified soldier. 
Soldiers are collectively grouped into "Forces" and act 
according to a user defined rule set for that force. That 
is, all agents of a Force are homogeneous with respect 
to their rule set definition. In order to create a truly 
heterogeneous battlefield of agents, each agent must be 
assigned to a Force consisting only of that agent. As this 
could result in the creation of a large numbers of Forces 
and is not generally practical, two opposing Forces are 
usually used. Additional to the rule sets, interaction in 
the form of weapons effects determines agent attrition. 
Agents are defined to exist in three health states, alive, 
injured and killed and can obviously only exist in one 
of these states. For each time step the position for each 
agent is determined according to the application of the 
previous time step data to the rule set. Following this 
positional update, attrition is calculated based on the 
defined hit/kill probabilities. Attrition is applied to a 
specified maximum number of enemy agents (from zero 
to all agents) within a specified range. 

Movement of ISAAC agents (ISAACAs) is determined 
by minimisation of a penalty function calculated at each 
time step. Proximity of other ISAACAs, both friendly 
and enemy, and proximity of each force's goals are de- 
fined as the six personality weights used to calculate the 
penalty at each location within the ISAACA's movement 
range rjv/. Note that the movement range and all other 
ranges used in ISAAC, is a square, not a circle, of radius 

Own force goals W5 are usually located at the starting 
side or corner of the domain/battlespace of their respec- 
tive force, with the opposition force goals wq located at 
the opposite side. As the proximity of an agent to the 
opposition goal increases, so the effect of the We term in- 
creases, increasing the speed at which it approaches the 
goal. Usually the value of W5 is set to zero, thus the own 
goal has no effect of an agent. 
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FIG. 1: Ranges used in ISAAC 



Personality weights: 
Wi alive friendly 
W2 alive enemy 
W3 injured friendly 
W4 injured friendly 



Ws friendly flag The form of the penalty function is: 

We enemy flag 
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where = 1,6 are personality weights as mentioned. 
The scaling factors V^?"/ and V^r^, number of ISAA- 
CAs Ni within sensor range and distances d[i; {x, y)] of 
those Ni ISAACAs form a discrete convolution, dnew and 
doid are the distances to the flags in question from each 
potential new position and from the original position re- 
spectively. 

At each time step this penalty is calculated for all po- 
tential moves an agent may make, including the penalty 
for remaining in the original position. When multiple 
new positions of equal minimum penalty occur, the new 
position is chosen randomly from this set. Interestingly, 
Ilachinski refers to these six personality weights as consti- 
tuting a local rule set, yet the influence of these weights, 
especially W5 and wg, can span the entire domain. 



A. Meta-Personality 

In addition to this penalty function, there are six ad- 
ditional rules that may be implemented, Advance Con- 
straint, Cluster Constraint, Combat Constraint, Mini- 
mum distance to friendly ISAACAs, Minimum distance 
to enemy ISAACAs, and Minimum distance to own flag. 
They are collectively termed a 'meta-personality' and 
modify the calculation of the penalty function. These are 
effectively variations on some of the personality weights 
of Eq [3] using user deflned threshold and constraint 
ranges. 

These meta-personalities will mean our PDEs as given 
in ([T]) and ^ will require modification in order to take 
these behaviours into account if we wish to make a di- 
rect comparison between our model and ISAAC. As these 
traits are not diffusion or reaction based, our fvei term 
will be the section modifled. 



1. Advance Constraint 

This constraint consists of specifying a threshold num- 
ber of friendly ISAACAs that must be within a given 
ISAACA's constraint range rc in order for that ISAACA 
to continue advancing toward the enemy flag. If the num- 



ber of friendly forces within the range exceed this thresh- 
old, the default weight +wq is used in Eq[3] If this is not 
the case, —wq is used so that the overall desire switches to 
a movement away from the enemy goal. Thus sufficient 
troop density is required in order for advancement. 



2. Cluster Constraint 

In order for the Advance constraint to be effective, 
a desire to be attracted toward friendly ISAACAs to 
form clusters is required. Again, the threshold number 
of ISAACAs required to be present within the constraint 
range rc is defined by the user. If this threshold is not 
met, the ISAACA will move in the direction of the high- 
est density of friendly forces calculated within the sensor 
range. Once this threshold is reached, an ISAACA will 
no longer move toward friendly ISAACAs, effectively set- 
ting the parameters wl = w3 — 0. 

Note that this constraint combined with wi and W3, 
mimics the attraction/repulsion kernels in our continuous 
equations. If the cluster constraint is not activated by the 
user, an artificial default lattice repulsion is present. Due 
to the density limitations of the lattice, only one agent 
of any type may occupy a lattice site at any time step, 
so that troop density cannot contract to a higher density 
than one troop per cell. In order to maintain an average 
inner density of less than one troop per cell, the cluster 
constraint must be used in conjunction with wi and ^3. 



3. Combat Constraint 

The Combat constraint is conceptually very similar to 
the Advance constraint, although this constraint is con- 
cerned with defining minimum conditions for engaging 
in combat with enemy ISAACAs. Two ranges are used 
to determine the number of friendly and enemy ISAA- 
CAs, N friendly (rc) and N^nemy{rs) respectively. For ad- 
vancement in the direction of greatest enemy concentra- 
tion within the range (rg), the threshold troop difference 
Ac = Nfriendiyirc) - NeneTay{rs) Hiust be cxcceded. If 
this threshold is not exceeded, movement is determined 
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using W2 = -W2,de fault and W4 = -w^^de fault, where 
W2,defauit and Wi^defauit are the default weights for mov- 
ing toward ahve and injured enemy IS A AC As. That 
is, when a numerical advantage is reached, the agents 
will advance along the gradient of highest enemy concen- 
tration, otherwise a numerical advantage has not been 
reached and the agents will retreat down the gradient of 
highest enemy concentration. 

Setting the combat threshold to a large positive num- 
ber gives a very defensive force, whereas for a large neg- 
ative value a force will pursue the enemy despite the nu- 
merical disadvantage. 

^. Other Meta-Personality and Agent Constraints 

Genetic algorithms, local command and other con- 
straints may be incorporated into ISAAC and are not 
discussed in this research. Published example scenarios 
either do not utilise these constraints or are unaffected 
upon their removal. 

III. COMPARISONS WITH ISAAC SCENARIOS 

During the prosecution of this research, one major 
shortcoming of using an agent based model of this na- 
ture was highlighted - stochasticity. Those results shown 
in the figures contained within the scenario descriptions 
of [6] are not necessarily indicative of any expected be- 
haviour. Although there is some facility for the collection 
of basic statistics such as average cluster size, spatial en- 
tropy or Red/Blue interpoint distance, there has been 
no undertaking to establish whether the given realisa- 
tions are true representatives_of their scenario parame- 
ters. Graham and Moyeed [12] note that Langrangian 
models are akin to experiments rather than any theo- 
retical undertaking and provide a framework for estab- 
lishing reliability of such results. We do not pursue the 
application of such a framework here. Rather we take 
the ISAAC results presented at face value while noting 
the frequency of realisations that match the stated be- 
haviour, and through our continuous ap pro ach provide a 
type of verification in a similar vein to [l3'|. All simula- 
tions use a domain of size 100 x 100 grid points. 

A. Classic Fronts Scenario 

1. Description 

We begin with the input file included with the Ein- 
stein Test Release Version 1.0.0.4 Beta, Build Date 2000, 
for the "einstein_classic_fronts" scenario which Ilachinski 
likens to a clash between two viscous fluids. The two 
loosely grouped forces collide and align in relatively sta- 
ble long thin fronts. Attrition between these fronts and 
the inherent randomness of the movement updates causes 



a discrepancy in density at either the upper or lower point 
of the formations. Once this occurs, the forces are able 
to slowly filter around each other and proceed to their 
respective goals. 



2. ISAAC Parameters and Results 

We now investigate the minimum number of param- 
eters and their values required to display the original 
behaviour. Firstly, the parameters for alive and in- 
jured friendly and enemy ISAACAs are made equal such 
that there is no distinction in behaviour of an alive or 
injured ISAACA. Personality weights towards friendly 
forces (and thus clustering effects) were set to zero and 
found not to affect the overall "front forming" behaviour. 
Switching the Advance parameter off also had no effect. 
The final parameter values used are in Table [D 
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FIG. 2: These snapshots were generated using the pro- 
vided input file and demonstrates equivalent behaviour 
to those in 



3. PDE Modification 

Our original equations are modified to reflect the addi- 
tion of the combat meta-personality. This is effectively a 
modification of the velocity term fvei as we wish for only 
the direction of movement to be affected by the presence 
of the opposite force. Diffusion, inter- force attraction and 
repulsion are to remain unaffected. 

fvel = V • {u(CuU + Aa{Ka * u) ~ ArU{Kr * u))} (4) 
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We propose the following form of the velocity term, tak- 
ing force u as an example; 

f C' if {JJj^ u dxdy - JJj^ v dxdy) > Ac 
C=<^ (5) 
I —C if ijjj^ u dxdy — JJj^ v dxdy) < Ac 

where (J udxdy) is the number of friendly forces and 
(J vdxdy) the number of enemy forces within the circu- 
lar domain with radius corresponding to the sensor range 
rg. C is a constant and represents the overall movement 
toward the u or friendly flag. If a numerical advantage 
greater than Ac is reached, the force will proceed with a 
constant velocity toward its own flag. Should this numer- 
ical advantage not be attained, the velocity of the force 
is reversed such that it will retreat away from its flag. 



4. PDE Results 
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FIG. 3: Our model: Classic Fronts. 



In order to introduce an asymmetry to the continuous 
model so as to mimic the asymmetry seen in the ISAAC 
results, the initial positions of each force distribution is 
changed such that the velocity vectors are offset. That 
is, the forces will no longer collide "head on" . 
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FIG. 4: Our Model: Classic Fronts comparison with 
forces initially offset. 
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ISAAC has two areas of inbuilt randomness - new po- 
sition selection when there are multiple positions with 
equal penalty function values, and the calculations of ca- 
sualties/fatalities. This randomness results in slight dif- 
ferences in the distribution of the two forces whereas the 
non-random continuous version has equal distributions. 
The ISAAC results in Figure [5] shows the forces "slip- 
ping" around each other due to these slight variations 
in distributions and then proceeding to their respective 
goals. Figure [3] shows that our continuous version also 
forms long thin fronts that remain stationary. Due to the 
steady reduction in both forces due to firing effects, their 
densities gradually decline however the position of the 
fronts remain stationary. These tactics represent a clas- 
sic style of attrition warfare. If the simulation is allowed 
to progress, both forces will eventually decline to zero. 
There is no randomness present to generate the slight dif- 
ferences or asymmetry in spatial distribution that leads 
to the forces manoeuvring around each other as in Fig- 
ure [H 



FIG. 5: Losses for Figure H] 

Again the formation of fronts is present however they 
form at a slight angle rather than vertically due to the 
initial offset in position. A slight oscillation in the move- 
ment of each force is visible in Figure [4] as they slowly 
rotate around each other. This is more noticeable in the 
movie of this scenario. Once the forces pass, circular 
profiles with the respective desired minimum densities 
are reformed. Comparing Figure [2] to Figure SI the tight 
troop formation of the ISAAC agents seen prior to the 
forces making contact in Figure [5] is not seen after the 
confrontation. In the ISAAC results, troops stream to 
their respective flag in single file rather than proceeding 
in formation due to the effect of increasing flag proxim- 
ity (wq) on the penalty function. However, maintaining 
a coherent profile throughout the entire simulation, es- 
pecially when a force is not in contact with another, is a 
highly desirable feature as a long thin profile may have 
an increased vulnerability to attack. 



B. Precess Scenario 

1. Description 

This scenario is described by Ilachinski as showing a 
simple example of an emergent behaviour. Red and blue 
forces are quite different in personahty, with red prefer- 
ring to remain in close proximity while blue actively seeks 
the red force within sensor range. Initially the rapidly 
advancing blue force is loosely formed and the slower 
red force advances in a higher density and thus smaller 
formation. As the forces enter into each other's sensor 
range, the blue force partially surrounds the red and a 
slow precession begins. The red force slowly continues 
advancing toward its flag, constantly pursued by the blue 
force as shown in Figure [H This precession behaviour, 
or rotation of both forces around an axis, as shown in 
the published ISAAC results and Figure [8] occurs in only 
a small percentage of the test runs. When precession 
was observed, both clockwise and anti-clockwise preces- 
sion occurred with the forces often colliding with the do- 
main boundaries. Offsetting the position of the blue flag 
slightly to the left or right did result in a higher frequency 
of precession observed, with precession directionality de- 
pended upon the flag offset. 

Without the flag offset the usual behaviour is as follows 
using the initial conditions as given in Table H As the 
forces come into sensor range, the blue force surrounds 
the red until the majority of the blue force is located on 
the south western side of the red force. As the person- 
ality weights for engaging in combat is stronger than for 
continuing to the blue flag, the blue force actively pur- 
sues the tightly clustered red force as it moves toward 
its red flag goal. Both forces remain in coherent forma- 
tions throughout the duration of the scenario as seen in 
Figure [6] 

We propose that, similar to the Classic Fronts scenario, 
this precession behaviour arises due to the spatial asym- 
metry of the forces. Depending on the spatial distribu- 
tion on falling into sensor range, the precession will be 
either clockwise or anticlockwise, and offsetting the ini- 
tial positions exacerbates this asymmetry thus increasing 
the frequency of precession occurring. 

2. ISAAC Parameters and Results 

Unlike the Classic Fronts scenario, no input file was 
provided by the developer of ISAAC with the software 
installation, only a series of snapshots with a partially 
complete list of parameter values. Those unspecified val- 
ues were determined through a trial and error procedure 
until the behaviour shown in Q could be reproduced. 
In a similar manner as for the Classic Fronts scenario, 
we set the Advance and Minimum Distance parameters 
to zero and find that precession behaviour is still pro- 
duced although more infrequently, suggesting that these 
parameters are unnecessary for this behaviour to occur. 
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FIG. 6: ISAAC Precess Scenario Screenshots 



Figure [6] shows the usual results found when running 
the input file as detailed in Table H] Note the absence 
of precession. As indicated above, precession in both 
directions can be observed with the direction determined 
by the spatial asymmetry generated in that particular 
instance. 





Fraction of Remaining ISAACAs 














/ Red (Total) 




1 




















Blue (Totai) 


















ft " 







FIG. 7: Losses for Figure [S] 
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FIG. 8: ISAAC Precess Scenario Screenshots with artifi- 
cial Offset 



By introducing an artificial offset, precession is almost 
always observed (Figure [S]) and is dictated by the direc- 
tion of the given offset. 



3. PDE Modifications 

We now seek to modify the velocity term in a similar 
way as with the Classic Fronts scenario. Our modification 
to the velocity term described in ([5]) does not adequately 
take into account the differences in offensive/defensive 
personalities as described in the ISAAC literature and 
must be expanded in order to include this. Firstly, we 
define the number of friendly {Ny) and enemy forces (iV„) 
within the circular domain of radius rs {Sensor range), 
used to determine the Combat constraint by: 



Ny = 



V dxdy; Ny — u dxdy; 



(6) 



In order to allow for the two different types of attacking 
personalities, a switch attack in the form of an integer of 
value 1 or —1 is included. This allows for a distinction 
between two force types - an aggressive and a defensive 
force. An aggressive force will move toward the enemy 
regardless of superiority in numbers, while a defensive 
force will advance only with superior numbers and will 
otherwise retreat. We now alter the velocity term again 
and arrive at the following form: 



C + JJj^v dxdy 



Nu- Ny> Ac 



C + attack x JJj^v dxdy ~ < Ac 



(7) 



With a numerical superiority above the given Ac, the 
force has an additional attraction up the gradient of 
greatest concentration of the opposing force. If this supe- 
riority level is not reached, this attraction is reversed, be- 
coming a repulsion down the gradient of greatest concen- 
tration. This effectively mimics the Combat constraint 
of ISAAC as described in Section [17X31 



4. PDE Results 

Figure [9] shows the first comparison at the ISAAC pre- 
cession results and is markedly similar to Figure [Sj Upon 
falling into the Sensor range of the smaller footprint 
higher density Red force, the larger footprint lower den- 
sity Blue force rapidly moves to surround it as attack = 1 
for Blue. As with the ISAAC scenario, this aggressive 
force pursues Red for the remainder of the simulation 
with the majority of the force located on the south west- 
ern side of the Red force and a small portion of Blue 
located on the north eastern side of the Red force. 
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FIG. 9: Our Model: Precess approximation. 
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FIG. 10: Losses for Figured] 

Note that no precession behaviour is seen which corre- 
sponds to the majority of ISAAC Precess scenario sim- 
ulation results. We now offset the initial positions of 
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the forces as we did with the ISAAC scenario such that 
they pass to the left or right of each other to ascertain 
whether both clockwise and anticlockwise precession be- 
haviour will be achieved in our continuous model. We 
begin by attempting to induce anticlockwise precession. 
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FIG. 11: Our Model: Precess comparison with initial 
positions offset. Anticlockwise precession. 
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FIG. 12: Losses for Figure [H] 

As expected, precession occurs due to the spatial asym- 
metry introduced into the scenario through the offset 
of the initial spatial distributions. The behaviour is 
markedly similar to that seen in Figure [H Upon falling 
into sensor range, the Red force is partially surrounded 
by Blue and an anti-clockwise precession occurs. The 
Red force is then closely pursued by the majority of the 
Blue force following and a small portion directly in the 
path of Red. Setting the initial positions with the reverse 
offset results in a clockwise precession. 

We propose that, similar to the Classic Fronts sce- 
nario, the precession behaviour seen in the ISAAC results 
(Figure [8|) arises due to the spatial asymmetry of the 
forces. Depending on the spatial distribution on falling 
into sensor range, the precession will be either clockwise 
or anticlockwise, and offsetting the initial positions of 
the forces exacerbates this asymmetry thus increasing 
the frequency of precession occurring. An artificial initial 



spatial asymmetry in the continuous version is necessary 
before precession is observed. 

As a result of the modification of the advection term 
through the addition of the attack variable, minimal 
overlap between the forces is observed which is in keeping 
with the ISAAC behaviour. However the positive value of 
attack for the Blue force is sufficient to give rise to pursu- 
ing of the Red force, yet the threshold value prevents any 
significant overlap. It is the effect of the repulsion down 
the gradient of highest concentration of enemy density 
which prevents any significant overlap. 



Circle Scenario 



Description 



This scenario shows strikingly similar behaviour to the 
Precess scenario. Both forces form dense formations with 
Blue slightly denser due to the higher cluster variable. 
The more aggressive Red force envelops the Blue force 
entirely upon coming into sensor range and it is at this 
point, time = 140, that the published snapshots stop. 
This can give the impression that this formation remains 
static or stable after this time. When running the simu- 
lation past this time step, the similarity with the Precess 
scenario from time — 120 onwards becomes apparent. 
Blue forces continue to move towards the Blue flag and 
Red forces begin to concentrate on the opposite side (Fig- 
ure [T31 time — 200). For the remainder of the simula- 
tion. Red forces pursue the Blue force toward the Blue 
flag. In this particular time series, a smaller section of 
the Blue force is separated and pursued by the Red force 
(time = 220) in the same manner as the larger section. 



2. ISAAC Parameters and Results 



Unlike the Classic Fronts scenario and similar to the 
Precession scenario, no input file was provided with the 
software installation, only a series of snapshots with a 
partial list of parameters and corresponding values. Ini- 
tial distribution position and size values were not pro- 
vided. Numerous test runs using the given parameters 
showed a less dense blue force advancing more rapidly 
than the higher density red force. Setting the Advance 
and Minimum Distance parameters to zero still produced 
the same behaviour. 
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FIG. 13: ISAAC Circle Scenario Screenshots 
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FIG. 15: Our Model: Circle approximation with differing 
fire parameters. 
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PDE Results 
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Using similar parameter values to the Precess scenar- 
ios and using the same form of the advection term ([7|. 
the following results are obtained. The only differencea 
between the two scenarios are the fire parameters. 
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FIG. 14: Our Model: Circle approximation. 



FIG. 16: Losses for Figure [Ml and Figure [T5l 

Similarly to the ISAAC results shown in Figure [T31 in 
FigurefTilthe denser Red force is surrounded by the more 
aggressive Blue force. While the repulsion due to force 
inferiority maintains a degree of separation between the 
two forces, some overlap is present at the edges. Once 
the Red force is surrounded, overall movement slows to 
almost a halt. Density is gradually lost due to attrition 
and there is a slight imbalance in the Blue force's distri- 
bution around Red. The slow advance of the Red force 
and imbalance in the Blue force distribution results in 
an eventual breakaway of Red. The higher attrition rate 
used in Figure [T5] shows that this almost stationary phase 
has a much shorter duration as attrition speeds the gener- 
ation of an imbalance in the surrounding force's density. 

Note the nonlinear form of density loss shown in Fig- 
ure [ini Lauren states that casualty rates are uneven on 
a turbulent nonlinear battlefield. This is held as an im- 
portant distinction between complex adaptive and con- 
vectional combat models. Here we have shown that a 
conventional model is indeed capable of producing inter- 
mittent density losses. 
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IV. DISCUSSION 

Cellular automata can be difficult to use for under- 
standing the underlying dynamics of combat as stochas- 
ticity can hinder the extraction of conclusions from a 
model. All scenarios presented here highlight the dan- 
gers associated with attributing intelligent reasoning to 
behaviour shown, when this can be explained quite sim- 
ply through the effects of the terms in our equations ^ 
and (21) and the spatial distribution of forces. This can be 
seen quite simply in the Classic Fronts scenario. Ilachin- 
ski describes the occurrence of the forces passing by one 
another as the agents having "found" a way to sneak 
around. We believe it is the differences in distribution 
caused the inherent randomness of ISAAC that results 
in agents located at the end of the distributions becom- 
ing sufficiently distanced so that repulsive forces from en- 
emy agents greatly reduces with respect to the attractive 
forces of the goal. This dominant attractive force then 
results in those agents progressing towards the goal. Al- 
though the attraction/repulsion of goals in our equations 
is not weighted with respect to separation distance as it 
is in the ISAAC penalty function, a slight offset or asym- 
metry of initial distributions demonstrates this same type 
of behaviour. 

This can also be seen clearly in the Precess scenario 
where the infrequently observed precession behaviour 
arises from asymmetries of the force profiles upon com- 
mencing combat, rather than being an unexplainable 
emergent behaviour. The majority of ISAAC Preces- 
sion simulations show behaviour similar to those obtained 
from the continuous equivalent without the initial dis- 
tribution offset as shown in Figure [H Again the ob- 
served asymmetries are due to the stochasticity of ISAAC 
present in the movement and attrition algorithms, and 
are responsible for the precession seen in Figure [5] and 
• When mimicked in our continuous equations through 
initially offset distributions, precession was observed. In 
this case the underlying behaviour without the effects of 
stochasticity highlighted the similarity between the Pre- 
cession and Circle scenarios. Again Ilachinski infers a 
degree of intelligence or emergence by stating that there 
are a few stray 'leakers ' and an occasional group of a few 
Blue IS A AC As that choosefs] to leave the main battle and 
head toward Red's flag. This is similar to the comments 
made for the Classic Fronts scenario and our explanation 
for that behaviour - differences in distribution result in 
agents becoming sufficiently distanced from enemy agents 
such that the goal terms in the penalty function become 
significant. 

Our deterministic approach encapsulates basic moti- 
vational factors and demonstrates a variety of spatial 
behaviours. Cohesive troop movement has only been 
achieved artificially in previous work: by using a desired 
initial distribution, low diffusion constant, and a suffi- 
ciently short overall simulation time, excessive diffusion 
is prevented and the troop profile cannot diffuse to a un- 
realistic spread. We have demonstrated that by using 



relatively simple and physically meaningful form of par- 
tial differential equations, cohesive troop movement can 
be achieved and maintained, even when suffering loss of 
density through fire. Forces and firing coefficients re- 
main homogeneous, a criticism of the traditional Lanch- 
ester approach, yet the nonlinear nature of the equa- 
tions are able to mimic those seen in ISAAC that have 
been labelled as complex. A continuum of forces is able 
to behave in a manner similar to a collection of indi- 
vidual autonomous agents, and shows decentralised self- 
organisation and adaptation of tactics to suit a variety of 
combat situations. This is a significant step toward de- 
veloping a set of realistic continuous equations for combat 
modelling. 

Lauren [§] states that complex adaptive models of com- 
bat, such as ISAAC or MAN A, display a rich variety 
of behaviour, a battlefield that is no longer linear and 
agent adaptivity. How exactly do the agents perform this 
adapting? The evolution of each agent's position is de- 
termined by a penalty function that remains unchanged 
throughout the entire simulation. There is a danger in 
the anthropomorphisation of agents, insinuating agents 
have a reasoning and planning ability when this is obvi- 
ously not the case. Also these types of wargames con- 
centrate heavily on the addition of extra communication 
ability between agents, shifting the emphasis to global or 
increasingly complex nonlocal features. MANA includes 
many more states and subsequently many more triggers 
need to be defined to facilitate switching between these 
states. Increasing the number of required parameters can 
cloud the process of deriving insight, a danger which has 
been shown in many of the basic scenarios presented here. 
For example the removal of parameters such as the Ad- 
vance constraint in the Precess scenario did not prevent 
the reproduction of similar precession behaviour. Lau- 
ren also states that conventional combat models behave 
largely as a series of attrition- driven fights and ISAAC 
entities will only fight if conditions are suitable. If our 
model is viewed as a conventional combat model as it 
is based on using Lanchester firing terms, it can be ar- 
gued that it also behaves as an ISAAC model due to the 
form of the Spatial Dynamics terms (fdiff and fvei)- Con- 
versely it could be argued that ISAAC is a type of con- 
ventional combat model as the movement of each agent is 
determined in a comparable way to our continuous model 
with the inclusion of randomness. 

Considering the popularity of cellular automata based 
wargames in military and complex adaptive systems re- 
search, it is imperative that complementary avenues of 
research are undertaken in order to gain a greater in- 
sight into the nature of combat. For a research topic 
such as combat modelling, emphasis of rare events or a 
misunderstanding of observed behaviour can have signifi- 
cant consequences. This research has demonstrated that 
seemingly emergent and intelligent behaviour of agents 
sometimes seen only in atypical scenarios, can be ex- 
plained through the numerical analysis of their contin- 
uous equivalent. The need for many multiple agent sim- 
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ulations and the application of data mining techniques 
becomes much reduced with the simuhaneous use of a 
continuous modeL MANA, much Uke ISAAC, can re- 
quire approximately 600 runs to establish a mean result 
for some scenarios [13]. Using the approach employed 
here it is much easier to establish a mean behaviour with 
our continuous form. 

By treating weapons systems in the simple forms as 
presented here, behaviour is less likely to be obscured by 
potentially highly nonlinear interaction terms or weapons 
effects. This could also lead to the exploration of the 
inherent nonlinearities in the Spatial Dynamics terms, 
which essentially describes the movement or tactics em- 
ployed by a force. For example, an increase in weapon 
lethality (kill probability) may not necessarily have the 
expected corresponding increase in casualties due to the 
effects of the enemy's tactics. It may not be sound to 
assume a doubling of this probability will yield double 
the number of casualties. 

We suggest that our continuous model be used in con- 
junction with agent-based wargames to act as a combined 
testbed for the purpose of concept exploration. 

V. FUTURE WORK 

This research did not address other ISAAC scenarios 
such as LOCALCMD or GLBALCMD where communi- 
cation between agents within a user-defined range affects 



the evolution of the scenarios. Inclusion of the local com- 
mand and global command functionalities and of com- 
mand personality in the penalty function ([3]) would re- 
quire further significant modification to our continuous 
model. Addition of these comparisons will result in a 
more complete development of a continuous counterpart 
to ISAAC. 

We have expanded our model to include a density re- 
sponse tactic where the equilibrium interior density is de- 
pendent upon the enemy density within the given sensor 
range. This is similar to predator avoidance modelling 
in CA fish schooling modelling in [l5| and has shown 
interesting dynamical changes in the test case scenarios 
thus far. We envisage publishing these results in the near 
future. 
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TABLE I: ISAAC Parameters for Scenario Screenshots. All other parameters are set to zero or No. 
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TABLE II: PDE Parameters. 
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